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Abstract 

Recently M. Ziman IjJ] criticized our approach for quantifyingthe required physical resources in the the- 
ory of Direct Characterization of Quantum Dynamics (DCQD) |2[] in comparison to other quantum process 
tomography (QPT) schemes. Here we argue that Ziman's comments regarding optimality, quantumness, 



Ti^j- ' and the novelty of DCQD are inaccurate. Specifically, we demonstrate that DCQD is optimal with respect 

o: 

\^ • to both the required number of experimental configurations and the number of possible outcomes over all 

"^ , known QPT schemes in the 2^" dimensional Hilbert space of n system and n ancilla qubits. Moreover, 



we show DCQD is more efficient than all known QPT schemes in the sense of overall required number 



j^ ■ of quantum operations. Furthermore, we argue that DCQD is a new method for characterizing quantum 

^: 

Q^. dynamics and cannot be considered merely as a subclass of previously known QPT schemes. 

> 



I. REPLY TO THE FIRST COMMENT: QUANTIFICATION OF RESOURCES 

In general, for a non-trace preserving completely positive (CP) quantum dynamical map acting 
on a rf dimensional quantum system, the number of independent elements to be characterized is 
exactly (i"^ ij]. In fact, in many important physical situations the quantum dynamical maps are 
effectively non-trace preserving. This phenomenon usually appears either as loss for photonic 
systems (due to inherent imperfections of optical elements), or leakage for atomic and spin-based 
quantum systems (due to interactions with photonic and/or phononic environments, spin-orbital 
coupling, etc.). Therefore, Ziman's statement, "Each quantum device acting on a d dimensional 
quantum system is described by d'^ — (P independent parameters that has to be specified in ar- 
bitrary (complete) process tomography scheme." yj] is not completely accurate. Even though 
one can always consider a larger Hilbert space such that the map becomes trace-preserving, this 
mathematical trick has little or no physical/practical significance, since in general we do not have 
full control over such an extended space. I.e., we cannot arbitrarily redefine our system and its 
environment in real life cases. 

In Ref. ^ we only compared our Direct Characterization of Quantum Dynamics (DCQD) 
scheme with separable quantum process tomography (QPT) methods, including standard quantum 
process tomography (SQPT) iS\ and the separable AAPT ^. Ziman raises the possibility of using 



non-separable QPT methods and proposed new resource quantification tables (the second and third 
tables in Ref. yJ]). As we argue below in a detailed analysis of such schemes, these tables are, 
unfortunately, incomplete and/or inaccurate. 

In general, the required state tomography in AAPT could also be realized by non-separable 
quantum measurements. These measurements can be performed either: (1) in the same Hilbert 
space, such as mutually unbiased bases (MUB) measurements O. la], or (2) in a larger Hilbert 
space, such as a generalized measurement or POVM I7J]. However, we find that these types of 
measurements would hardly have any practical relevance in the context of QPT, because they 
require many-body interactions that are not experimentally available. In the next two subsections 
we address each of these two approaches separately. 



A. AAPT with Mutual Unbiased Bases Measurements 

The AAPT scheme utilizes the degrees of freedom of an auxiliary system, B, in order to char- 
acterize an unknown quantum dynamical map acting on a principle system A. The information 
about the dynamics is obtained by complete quantum state tomography of the combined system 
and ancilla. Quantum state tomography in itself is the task of characterizing an unknown quantum 
state by measuring the expectation values of a set of non-commuting observables on the subensem- 
ble of quantum systems prepared in the same state. For characterizing the density operator of a 
d-dimensional quantum system, there are in general (P — 1 non-commuting observables to be 
measured. The minimal number of non-commuting measurements, that corresponds to a mutually 
unbiased basis, is d + 1 (for systems with d being a prime or a power of prime) O. |6|]. A set of 
bases, in a given Hilbert space, are mutually unbiased if the inner products of each pair of elements 
in these bases have the same magnitude. 

Let us consider the case of characterizing a non-trace preserving dynamical map acting on a 
single qubit A, using a single ancilla qubit B. For such a two-qubit system d = A, and the number 
of MUB for the required state tomography is five. Therefore, the minimum number of ensemble 
measurements (experimental configurations) in the AAPT scheme in a single qubit case is five, 
(as opposed io four mentioned in Ref. [1]). The first measurement provides four independent 
outcomes and the last four measurements each provide three independent outcomes. I.e., we have, 
1 X [(1 X 4) + (4 X 3)] = 16, where in each term in the sum the first number (1,4) represents 
the number of required measurements per input state, and the second number (4, 3) represents the 
number of outcomes for each measurement; 16 corresponds to the total number of independent 
outcomes. This should be compared with 4 x (1 x 4) = 16 in the DCQD scheme ||2|]. The first 
number (4) is the number of required different input states (compared to 1 in the AAPT scheme). 

It should be noted that if we know the local state of the ancilla (i.e., if we know the results of 
the measurements I^ ® I^ , I"^ ® a^ , I^ ® Uy and /^ ® erf from our prior knowledge about the 
preparation and the fact that the output state has trace unity), then we need to find only d'^ — d^ 
parameters for the superoperator. But even in this case, we still need^ve different measurement 
setups. For a simple proof, let us examine the MUB of a two-qubit system, see the table in Fig.l 
(also Ref. ^). Obviously, for a non-trace preserving map we need at least five (ensemble) mea- 
surements corresponding to each row (measuring simultaneously the commuting operators in the 
first two columns). Now the question is, how many measurements are needed if we already know 
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FIG. 1: A mutually unbiased basis (MUB) of the 4-dimensional Hilbert space of the principal qubit (A) and 
its ancilla (B). Each row represents a different experimental configuration. 

the local state of the ancilla? The answer becomes clear if we note that the measurements I^0a^, 
I"^ ® (Jy, and I"^ ® af appear in the second column of the first three rows and are redundant. How- 
ever, we still need to perform all three (ensemble) measurements corresponding to the first three 
rows, since the measurements in the first column (corresponding to the local state of the principal 



qubit \- a^®I ,ay®I and ay ® / ) do not commute. These three measurements plus the two 
measurements related to the 4th and 5th rows (corresponding to measuring the correlations of the 
principal qubit and the ancilla) result in five measurements overall. Note that the above argument 
is independent of the basis chosen, since in any other basis the measurements corresponding to the 
local state of the ancilla always appear in different rows due to the non-commuting properties of 
the Pauli operators. 

For the case of n qubits, and by using n ancillary qubits, the required measurements for the 
AAPT scheme can be performed using three different types of methods: (a) using 16" (separable) 
joint single-qubit measurements on the n principal and n ancillary qubits, or (b) using 5" mutually 
unbiased bases measurements (tensor product of MUB measurements on two-qubit systems), or 
more efficiently (c) using 4" + 1 (non-separable) mutually unbiased bases measurements on the 
Hilbert space of all 2n qubits. The latter method requires many-body interactions between all 
2n qubits, which are obviously not naturally available. The required global Hamiltonian can be 

2n 2n 2n 

written as // = X) Oj = J] (^ (^au ky where [Oj, Oj] = for any i,j G {1,2,..., 2n}, such 



j=i 



j=i fc=i 



that their common eigenvectors form a MUB, and cr^^-j,-, is a Pauli operator (a^, ay or a^) acting 
on the k^^ qubit. I.e. one should simultaneously measures 2n commuting Hermitian operators 



Oi E {01,021 •••02n}; each operator 0/ in itself acts on 2n physica' 
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qubits as Oi={ai 



"2,03,. ■■} 



'^ai'^a2^a3---'^a2n (^"S-' f°^ ^^c casc of thrcc qubits see Fig. 2 of Ref. ||6|]). 

The operators Oj and Oj commute globally and are made of tensor products of Pauli operators, 



however they cannot be simuhaneously measured locally, i.e. by using only single-qubit devices. 
The reason is that according to the Heisenberg uncertainty principle, the outcome of each local 
measurement a^ for the operator Oi completely destroys the outcome of measuring u^ ^ a\ for 
other operators Oj . In principle, one could simulate the required many -body interactions in AAPT 
(for each of the 4" + 1 measurements) by a quantum circuit comprising about O {v?) single and two- 
qubit quantum operations (with the assumption of realizability of non-local two-body interactions, 
i.e., with having access to two-qubit interactions between every pairs of the 2n qubit system) 
||8|]. For a simple proof, we note that the measurement of an operator in the form a].cr1a^^...a1^ 
requires 2n sequential CNOT operations. For measuring a more general operator of the form 
Oj={Qi, 02,03 ..} = ^ai'^Q2'^Q3---^a2n' ^^ '^^^^ ^'^ additional 0{n) local single qubit operations to 
make an appropriate change of basis. Therefore, for measuring 2n operators Oj, one at least needs 
to realize (2ri)^ quantum operations. However, if only local two-body interactions are available 
(i.e., if we are restricted to using nearest neighbor interactions) then O(n^) single- and two-body 
quantum operations would be required. The reason is that the overall number of operations grows 
by a factor of 0{n) due to the cost of transporting each two-qubit gate. A modified table for 
comparing physical resources in different QPT methods is presented in Fig. 2. 

Note that the AAPT plus MUB measurements with 0{n^) [or 0{n^)] two-body interactions 
is essentially in the same complexity class as a quantum Fourier transform. This should be com- 
pared to the DCQD scheme with a single step of n CNOT (between each qubit and its ancilla) 
and n Hadamard operations to realize a single Bell-state measurement. In the context of estimat- 
ing quantum dynamics, the implementation of 0{n?) or O(n^) gates in AAPT is inefficient, since 
these operations not only increase the time execution of each measurement, they also create addi- 
tional errors that would be very difficult to discriminate from the actual quantum dynamical map. 
Moreover, the overall number of repetitions for these operations scales poorly with a given desired 
precision. 

Although MUB is not the most general measurement for state tomography of a 2ri-qubit sys- 
tem, it is well-understood (for systems whose dimensions is a power of prime) to be the optimal 
measurement scheme for such a task. Therefore, any other measurement strategy within the same 
Hilbert space results in more ensemble measurements than MUB. For other systems (whose di- 
mensions is not a power of prime) the scaling of the AAPT measurements becomes even worse; 
since in those systems, the existence and construction of MUB is not fully understood and there- 
fore in general one has to measure a complete operator basis of 2n-qubit system which has 2''" — 1 
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FIG. 2: Comparison of the required physical resources for characterizing an arbitrary non-trace preserving 
CP quantum dynamical map on n qubits. The overall number of experimental configurations (5th column) 
is obtained from multiplying the number of required input states and the number of non-commuting mea- 
surements at the output states. In principle, the required many-body interactions for AAPT-i-MUB can be 
simulated by a quantum circuit comprising about 0(?i^) [0{n^)] single and two-qubit quantum operations 
with the assumption of non-local [local] two-body interactions. However, the many-body interactions for 
a POVM in this context cannot be efficiently simulated. The overall number of quantum operations (9th 
column) is defined as the product of the number of experimental configurations and the number of single- 
and/or two-body operations required for simulating each measurement. Finally, the overall number of quan- 
tum operations for a desired precision is obtained by multiplying by a factor of 2^^~^''^'N the 9th column, 
where 2^"" is the number of possible outcomes in each QPT strategy, and N is the number of repeated mea- 



surements to obtain the precision e 



We note that DCQD is more efficient than all known QPT 



schemes with respect to the overall number of quantum operations (see 9th and 10th columns). 



members. 

In principle, one could devise intermediate strategies for AAPT, using different combinations 
of single-, two-, and many-body measurements. The number of measurements in such methods 
ranges from 4" + 1 to 16", which is always larger than what is required in DCQD ||2|], using 4" 
Bell-state measurements. Therefore, in the 2^" dimensional Hilbert space of the 2n system and 
ancillary qubits, DCQD requires fewer experimental configurations than all other QPT schemes. 



Using DCQD one can in principle transfer log2 2^" bits of classical information between two 
parties, Alice and Bob, which is optimal according to the Holevo bound ^. Alice can realize 
this task by encoding a string of 2^" bits of classical information into a(n) (engineered) quantum 
dynamics (e.g., by applying one of 2^" unitary operator basis to the n qubits in her possession 
and then send them to Bob). Bob can decode the message by a single measurement on 2n qubits 
using DQQD scheme M\. I.e., the overall number of possible independent outcomes in each 
measurement in DCQD is 2^", which is exactly equal to the number of independent degrees of 
freedoms for a 2n qubit system, therefore, a maximum amount of information can be extracted in 
each measurement in DCQD, which cannot be improved by any other possible QPT strategies in 
the same Hilbert space. 

B. A APT with Generalized Measurements 

In principle, it is possible to perform the required quantum state tomography at the output 
states of an AAPT scheme by utilizing a single POVM or generalized measurement M\. How- 
ever, for characterizing the dynamics on n qubits the number of required ancillary qubits should 

n 

be increased from n to 3n (see Fig. 3. of Ref |7]). This can be easily understood according 
to the Holevo bound. For extracting complete information about a quantum dynamical map (en- 
coded by 2*'^ independent parameters of the superoperator) in a single measurement, one needs 
a Hilbert space of dimension at least 2"^"; otherwise the Holevo bound cannot be satisfied. There 
are two major disadvantages of using such a POVM compared to all other QPT schemes. (1) The 
POVM measurement requires a general many-body interaction between 2n qubits that cannot be 
efficiently simulated. I.e., it requires an exponential number of single- and two-qubits quantum op- 
erations. (2) The number of required repetition of each measurement to obtain a desired precision 
grows exponentially with n. 

According to the general setting in Ref. IlZD, in order to implement a POVM for extracting all 
the information about any observable of an n-qubit system, one needs to realize a global normal 
operator H (a single universal quantum observable) in the Hilbert-Schmidt space of the principle 
system and an ancilla system in the form of H = J2'i=i Ei ® Pi. Here [H, H^^] = 0, {Ei] is an 
operator basis for the n-qubit Hilbert space of the system, and [Pi] is a set of projections over 
the ancilla Hilbert space (e.g.. Pi = \i){i\, where {\i)} is an orthonormal basis in the ancillary 
Hilbert space). The operator H has the most general form of an operator-Schmidt decomposition 



m 



llOll and cannot be simulated in a polynomial number of steps. It is known that in general at 
least 0(4^") single- and two- qubi t operations are needed to simulate such general many-body 
operations acting on 2n qubits llllll (see also llioll for different measures of complexity of a given 
quantum dynamics). 

One important disadvantage of all QPT schemes in a larger Hilbert space (that rely on calcu- 
lating 2'^" different joint probability distributions) is that each measurement has to be repeated 
by a factor 2^'^"^^'^ in order to build the same statistics as a SQPT scheme. Let us consider an 
implementation of SQPT with a desired precision of e ~ -4^ in characterizing parameters of a 
superoperator, where e represents the standard deviation and N is the number of repeated mea- 
surements. Since each measurement in SQPT has 2" possible outcomes, the precision e can be 
obtained by A^ = ^ measurements. Note that in order to obtain a similar statistical error e with 
other methods of QPT, with 2'^" possible outcomes, we need to perform A^' = 2^''^^^'"'N measure- 
ments for each experimental configuration. Therefore, the actual number of measurements for a 
POVM strategy, with 2"^" possible outcomes, grows by a factor of 2^" with respect to SQPT and 
2^" with respect to DCQD - see the last column of the table in Fig. 2. We note that the overall 
number of quantum operations is still optimal for DCQD for any desired precision. 

n. REPLY TO THE SECOND COMMENT: INDEPENDENCE FROM AAPT AND USAGE OF 
QUANTUMNESS 

Here, we argue that the DCQD scheme is an independent algorithm from the AAPT scheme. 
First, we note that DCQD has a different scaling from the AAPT scheme in the sense of the overall 
number of experimental configurations, or similar scaling while using only two-body interactions 
compared to many-body interactions in the non-separable AAPT. Clearly, such results could not 
have been obtained if DCQD were merely a subclass of AAPT. 

Second, the required input states for DCQD must be entangled, which is in complete contrast 
to AAPT, where entanglement is not required at the input level. E.g., for characterizing quantum 
dynamical populations ||2|], the input state in DCQD must be maximally entangled in order to 
form a nondegenerate stabilizer state (with two independent stabilizer generators). This follows 
from the quantum Hamming bound, according to which only nondegenerate stabilizer states can 
be utilized for obtaining full information about the nature of all 2^ different error operator basis 
elements, acting on a single qubit of a two-qubit system [jj]. For characterizing quantum dynamical 
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coherence, the entanglement in input states is absolutely necessary, for otherwise the expectation 
values of the normalizers always vanish and therefore do not provide any information about the 
dynamics. In addition, the error-detection measurements in DCQD (e.g., Z^Z^ , X^Z^) are 
fundamentally non-separable, which is again in contrast to AAPT which also can be performed by 
joint single-qubit measurements. 

Third, the DCQD method utilizes a different methodological approach to quantum dynamical 
characterization than the AAPT schemes. DCQD utilizes a set of entangled input states and pro- 
vides the set of commuting observables to maximize the amount of classical information about the 
dynamics, log2 d? ^ that can be obtained at each output state (according to the Holevo bound), with- 
out completely characterizing any of the output states. However, AAPT utilizes a single faithful 
(not necessary entangled) input state and provides the minimal set of non-commuting observables 
that should be measured in order to completely characterize the output state of the combined sys- 
tem and ancilla, and the amount of classical information that can be obtain in each measurement, 
\og2{d^ — 1), is less than the maximum allowable according to the Holevo bound (log2 d^)- There- 
fore, there cannot be any direct correspondence between DCQD and AAPT in any fixed Hilbert 
space of the system and ancilla. We believe that the only true similarity between AAPT and DCQD 
is the fact that both methods utilize the degrees of freedoms of an auxiliary system. 

III. REPLY TO THE THIRD COMMENT: DIRECT CHARACTERIZATION OF DYNAMICS 

In order to remove any ambiguity, we would like to define what we mean by direct charac- 
terization of quantum dynamics. We define a QPT method to be a direct method if it satisfies 
these two conditions: (1) It should not rely on complete state tomography of the output states. 
(2) Each experimental outcome (joint probability distribution of observables) give direct informa- 
tion about either a single element of the superoperator, {Xmm) or a specific known subset of the 
superoperator's elements (e.g., Xmm, Xnn,Re(xmn))- 

AAPT is not a direct method because it does not satisfy the condition (1). One could argue that 
when the local state of the ancilla is known {d^ parameters), only d^ — d^ additional parameters 
must be characterized. However, since in this case one eventually also has access to d"^ parameters 
{d^ of which were known from the beginning, the rest measured), this should clearly also count as 
complete state tomography 



IV. CONCLUSION 

We agree that demonstrating an absolute advantage of one QPT scheme over the other QPT 
methods requires a complete quantification of the complexity of the preparations and measure- 
ments procedures. Indeed, we present a more detailed analysis in a separate publication L^. In 
conclusion, we believe the following statements are true: 

1. DCQD is a quantum algorithm for complete and direct characterization of quantum dynam- 
ics, which does not require state tomography. 

2. We have proved that DCQD is optimal in the sense of both the required number of exper- 
imental configurations and the number of possible outcomes, over all other known QPT 
schemes in a given Hilbert space. 

3. DCQD is quadratically more efficient than all separable QPT schemes in the number of 
experimental configurations. 

4. A similar scale-up in the number of experimental configurations is achievable with the 
AAPT scheme and MUB measurements only if many-body interactions are realized (or 
simulated with 0{'n?) or O(n^) single- and two-body gates). 

5. In principle, by utilizing POVMs, a single experimental configuration is sufficient for a com- 
plete QPT, however, one should realize many-body interactions that are not experimentally 
available and cannot be efficiently simulated by single- and two-body interactions. More- 
over, the POVM strategy has to be repeated as many as 4" times more than DCQD to obtain 
a similar precision. 

6. DCQD is new method for QPT, and cannot be considered merely as a subclass of any known 
QPT methods. 

7. DCQD is the first theory that utilizes quantum error detection methods in quantum process 
tomography. 

V. APPLICATIONS AND FUTURE WORK 

We believe that a potentially important advantage of DCQD is for use in partial characterization 
of quantum dynamics, where we cannot afford or do not need to carry out a full characterization of 
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the quantum system under study, or when we have some a priori knowledge about the dynamics. 
We have already presented two examples in connection with simultaneous measurement of Ti and 
T2, and realization of generalized quantum dense coding tasks. Other implications and applica- 
tions of DCQD remain to be investigated and explored, specifically, for obtaining a polynomial 
scale-up in physical resources for partial characterization of quantum dynamics. We believe that 
in some specific regimes DCQD could have near-term applications (within the next 5-10 years) 
for complete verification of small quantum information processing units (fewer than five qubits or 
so), especially in trapped-ion and liquid-state NMR systems. For example, the number of required 
experimental configurations for systems of 3 or 4 physical qubits is reduced from 5000 and 65000 
(in SQPT) to 64 and 256, respectively, in our scheme. Complete characterization of such dynam- 
ics would be essential for verification of quantum key distribution procedures, teleportation units 
(in both quantum communication and scalable quantum computation), quantum repeaters, quan- 
tum error correction procedures, and more generally, in any situation in quantum physics where 
a few qudits have a common local bath and interact with each other. Another interesting exten- 
sion of DCQD is to develop a theory for closed-loop and continuous characterization of quantum 
dynamics by utilizing weak measurements for our error-detection procedures. 
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